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NLO QCD corrections to WW+je-t production at hadron colliders 
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We report on the calculation of the next-to-leading order QCD corrections to the production of 
W-boson pairs in association with a hard jet at the Tevatron and the LHC, which is an important 
source of background for Higgs and new-physics searches. The corrections stabilize the leading-order 
prediction for the cross section considerably, in particular if a veto against the emission of a second 
hard jet is applied. 
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INTRODUCTION 

The search for new-physics particles — including the 
Standard Model Higgs boson — will be the primary task 
in high-energy physics after the start of the LHC that is 
planned for 2008. The extremely complicated hadron col- 
lider environment does not only require sufficiently pre- 
cise predictions for new-physics signals, but also for many 
complicated background reactions that cannot entirely be 
measured from data. Among such background processes, 
several involve three, four, or even more particles in the 
final state, rendering the necessary next-to-leading-order 
(NLO) calculations in QCD very complicated. This prob- 
lem lead to the creation of an "experimenters' wishlist for 
NLO calculations" P, Q that are still missing for success- 
ful LHC analyses. The process pp — > W + W~+jet+X 
made it to the top of this list. 

The process of WW+jet production is an important 
source for background to the production of a Higgs boson 
that subsequently decays into a W-boson pair, where ad- 
ditional jet activity might arise from the production or a 
hadronically decaying W boson. WW+jet production de- 
livers also potential background to new-physics searches, 
such as supersymmetric particles, because of leptons and 
missing transverse momentum from the W decays. Last 
but not least the process is interesting in its own right, 
since W-pair production processes enable a direct precise 
analysis of the non-abelian gauge-boson self- interactions, 
and a large fraction of W pairs will show up with addi- 
tional jet activity at the LHC. 

In this letter we report on the first calculation of the 
process pp -> W+W^+jet+AT in NLO QCD. 
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FIG. 1: Two representative LO diagrams for the partonic 
process uu — > W + W~g. 




FIG. 2: Pentagon diagrams for the partonic process uu 

w+w-g. 



nels are related by crossing symmetry to the amplitude 
— > W + W - <7<7g. Two representative LO diagrams for 
the process uu — > W + W~g are shown in Figure [T] 

In order to prove the correctness of our results we have 
evaluated each ingredient twice using independent calcu- 
lations based — as far as possible — on different methods, 
yielding results in mutual agreement. 



Virtual corrections 



DETAILS OF THE NLO CALCULATION 

At leading order (LO), hadronic WW+jet produc- 
tion receives contributions from the partonic processes 
qq -► W+W-g, qg -► W + W~ 9 , and qg -> W+W~g, 
where q stands for up- or down-type quarks. Note 
that the amplitudes for q = u, d are not the same, 
even for vanishing light quark masses. All three chan- 



The virtual corrections modify the partonic processes 
that are already present at LO. At NLO these corrections 
are induced by self-energy, vertex, box (4-point), and 
pentagon (5-point) corrections. For illustration the pen- 
tagon graphs, which are the most complicated diagrams, 
are shown in Figure [2]for a partonic channel. At one loop 
WW+jet production also serves as an off-shell continua- 
tion of the loop-induced process of Higgs+jet production 
with the Higgs boson decaying into a W-boson pair. In 
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FIG. 3: Some diagrams with closed quark loops for the par- 
tonic process uu — > W + W~g. 



this subprocess the off-shell Higgs boson is coupled via 
a heavy-quark loop to two gluons; a sample graph for 
this mechanism is shown in Figure [3] together with some 
other typical graphs with a closed quark loop. 

Version 1 of the virtual corrections is essentially ob- 
tained as for the related processes ofttH andtt+jet Q 
production. Feynman diagrams and amplitudes are gen- 
erated with FeynArts 1.0 [5] and further processed with 
in-house Mathematica routines, which automatically cre- 
ate an output in Fortran. The IR (soft and collinear) sin- 
gularities are treated in dimensional regularization and 
analytically separated from the finite remainder as de- 
scribed in Refs. 0, @. The pentagon tensor integrals are 
directly reduced to box integrals following Ref. [7] . This 
method does not introduce inverse Gram determinants in 
this step, thereby avoiding numerical instabilities in re- 
gions where these determinants become small. Box and 
lower-point integrals are reduced a la Passarino-Veltman 
Q to scalar integrals, which are either calculated analyt- 
ically or using the results of Refs. BE EH- Sufficient 
numerical stability is already achieved in this way, but 
further improvements with the methods of Ref. [12j | are 
in progress. 

Version 2 of the evaluation of loop diagrams starts 
with the generation of diagrams and amplitudes via Feyn- 
Arts 3.2 [l3| which are then further manipulated with 
FormCalc 5.2 [14J and eventually automatically trans- 
lated into Fortran code. The whole reduction of tensor 
to scalar integrals is done with the help of the LoopTools 
141 ] . which also employs the method of Ref. [7J 



library 

for the 5-point tensor integrals, Passarino-Veltman |8| 
reduction for the lower-point tensors, and the FF pack- 
age 0,(3 

for the evaluation of regular scalar integrals. 
The dimensionally regularized soft or collinear singular 
3- and 4-point integrals had to be added to this library. 
To this end, the explicit results of Ref. @ for the vertex 
and of Ref. [13] for the box integrals (with appropriate 
analytical continuations) are taken. 



Real corrections 

The matrix elements for the real corrections are given 
by — > W + W _ gggg and — * W + W~ qqq'q' with a large 
variety of flavour insertions for the light quarks q and 
q' . The partonic processes are obtained from these ma- 
trix elements by all possible crossings of quarks and glu- 
ons into the initial state. The evaluation of the real- 



emission amplitudes is performed in two independent 
ways. Both evaluations employ (independent implemen- 
tations of) the dipole subtraction formalism [l8[ for the 
extraction of IR singularities and for their combination 
with the virtual corrections. 

Version 1 employs the Weyl-van-der-Waerden for- 
malism (as described in Ref. [19() for the calculation 
of the helicity amplitudes. The phase-space integration 
is performed by a multi-channel Monte Carlo integra- 
tor [13] with weight optimization 21 1 written in C++, 
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which is constructed similar to RacoonWW 
The results for cross sections with two resolved hard 
jets have been checked against results obtained with 
Whizard 1.50 Q and Sherpa 1.0.8 (25|. Details on 
this part of the calculation can be found in Ref. [2q ]. 
In order to improve the integration, additional channels 
are included for the integration of the difference of the 
real-emission matrix elements and the subtraction terms. 

Version 2 is based on scattering amplitudes calcu- 
lated with Madgraph [27J generated code. The code has 
been modified to allow for a non-diagonal quark mixing 
matrix and the extraction of the required colour and spin 
structure. The latter enter the evaluation of the dipolcs 
in the Catani-Seymour subtraction method. The evalua- 
tion of the individual dipoles was performed using a C++ 
library developed during the calculation of the NLO cor- 
rections for tt+jet [H. For the phase-space integration a 
simple mapping has been used where the phase space is 
generated from a sequential splitting. 



NUMERICAL RESULTS 

We consistently use the CTEQ6 [H \%k set of P ar " 
ton distribution functions (PDFs), i.e. we take CTEQ6L1 
PDFs with a 1-loop running a s in LO and CTEQ6M 
PDFs with a 2-loop running a s in NLO. We do not in- 
clude bottom quarks in the initial or final states, because 
the bottom PDF is suppressed w.r.t. to the others; outgo- 
ing bb pairs add little to the cross section and can be ex- 
perimentally further excluded by anti-b-tagging. Quark 
mixing between the first two generations is introduced 
via a Cabibbo angle 9c = 0.227. In the strong coupling 
constant the number of active flavours is N-p = 5, and 
the respective QCD parameters are Ag ° = 165 MeV and 
A^ s = 226 MeV. The top-quark loop in the gluon self- 
energy is subtracted at zero momentum. The running 
of a s is, thus, generated solely by the contributions of 
the light quark and gluon loops. The top-quark mass is 
m t = 174.3 GeV, the masses of all other quarks are ne- 
glected. The weak boson masses are Mw = 80.425 GeV, 
M z = 91.1876 GeV, and M H = 150 GeV. The weak 
mixing angle is set to its on-shell value, i.e. fixed by 
= 1 — s 2 ^ = Af^/M|, and the electromagnetic cou- 
pling constant a is derived from Fermi's constant — 
1.16637 x 10- 5 GeV" 2 according to a = V^G^M^sl/ir. 
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pp -> W+W~+jet+X 
Vi=1.96TeV 



LO (CTEQ6L1) 
NLO (CTEQ6M) 




PTjet > 20GeV- 
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FIG. 4: Scale dependence of the LO and NLO cross sections 
for WW+jet production at the Tevatron, where the renor- 
malization and factorization scales are set equal to p. 



We apply the jet algorithm of Ref. [30( with R = 1 
for the definition of the tagged hard jet and re- 
strict the transverse momentum of the hardest jet by 
PTjct > PT,jot,cut. In contrast to the real corrections 
the LO prediction and the virtual corrections are not 
influenced by the jet algorithm. In our default setup, 
a possible second hard jet (originating from the real 
corrections) does not affect the event selection, but 
alternatively we also consider mere WW+jet events with 
"no 2 nd separable jet" where only the first hard jet is 
allowed to pass the prjet cut but not the second. 

Figures [4] and [5] show the scale dependence of the inte- 
grated LO and NLO cross sections at the Tevatron and 
the LHC, respectively. The renormalization and factor- 
ization scales are identified here (/i = /i rC n = Mfact), and 
the variation ranges from \i = 0.1 Mw to p = 10 Mw- 
The dependence is rather large in LO, illustrating the 
well-known fact that the LO predictions can only provide 
a rough estimate. At the Tevatron the qq channels dom- 
inate the total pp cross section by about 90%, followed 
by the qg and qg channels with about 5% each. Scaling 
the renormalization and factorization scales simultane- 
ously by a factor of 4 (10) changes the cross section by 
about 70% (100%). At the LHC, the qg channels com- 
prise about 56%, followed by qq with about 28%. Surpris- 
ingly the scale dependence is much smaller than at the 
Tevatron: varying the scales simultaneously by a factor 
of 4 (10) changes the cross section by about 25% (50%). 

At the Tevatron (Figure [4]), the NLO corrections signif- 
icantly reduce the scale dependence for PTjct > 20 GeV 
and 50 GeV. We observe that around fx ~ Mw the NLO 
corrections are of moderate size for the chosen setup. At 
the LHC (Figure , only a modest reduction of the scale 
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pp -> W+W~+jet+X 
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FIG. 5: Scale dependence of the LO and NLO cross sections 
for WW+jet production at the LHC, where the renormaliza- 
tion and factorization scales are set equal to /i. 



dependence is observed in the transition from LO to NLO 
if W pairs in association with two hard jets are taken into 
account. This large residual scale dependence in NLO, 
which is mainly due to gg-scattering channels, can be 
significantly suppressed upon applying the veto of hav- 
ing "no 2 nd separable jet" . The contribution of the gen- 
uine WW+2jets events, which represents the difference 
between the two NLO curves in the plots of Figure 03 is 
also reduced if the cut on prjet is increased from 50 GeV 
to 100 GeV. The relevance of a jet veto in order to sup- 
press the scale dependence at NLO was also realized [31| 
for genuine W-pair production at hadron colliders. 

Finally, we show the integrated LO and NLO cross sec- 
tions as functions of px,jet,cut i n Figure [5] The widths of 
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FIG. 6: LO and NLO cross sections for WW+jet production 
at the Tevatron and LHC as function of Pt jet, cut- 

the bands, which correspond to scale variations within 
Mw/2 < [i < 2Afw, reflect the behaviour discussed 
above for fixed values of pt jot, cut- For Tevatron the re- 
duction of the scale uncertainty is considerable, for the 
LHC it is only mild unless WW+2jets events are vetoed. 
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